Zn2+ and Cu2+ Interaction with the Recognition Interface of ACE2 for SARS-CoV-2 Spike Protein.

The spike protein (S) of SARS-CoV-2 is able to bind to the human angiotensin-converting enzyme 2 (ACE2) receptor with a much higher affinity compared to other coronaviruses. The binding interface between the ACE2 receptor and the spike protein plays a critical role in the entry mechanism of the SARS-CoV-2 virus. There are specific amino acids involved in the interaction between the S protein and the ACE2 receptor. This specificity is critical for the virus to establish a systemic infection and cause COVID-19 disease. In the ACE2 receptor, the largest number of amino acids playing a crucial role in the mechanism of interaction and recognition with the S protein is located in the C-terminal part, which represents the main binding region between ACE2 and S. This fragment is abundant in coordination residues such as aspartates, glutamates, and histidine that could be targeted by metal ions. Zn2+ ions bind to the ACE2 receptor in its catalytic site and modulate its activity, but it could also contribute to the structural stability of the entire protein. The ability of the human ACE2 receptor to coordinate metal ions, such as Zn2+, in the same region where it binds to the S protein could have a crucial impact on the mechanism of recognition and interaction of ACE2-S, with consequences on their binding affinity that deserve to be investigated. To test this possibility, this study aims to characterize the coordination ability of Zn2+, and also Cu2+ for comparison, with selected peptide models of the ACE2 binding interface using spectroscopic and potentiometric techniques.

Allosteric communication between ACE2 active site and binding interface with SARS-CoV-2.

SARS-CoV-2, the virus causing COVID-19, initiates cell invasion by deploying a receptor binding domain (RBD) to recognize the host transmembrane peptidase angiotensin-converting enzyme 2 (ACE2). Numerous experimental and theoretical studies have adopted high-throughput and structure-guided approaches to (i) understand how the RBD recognizes ACE2, (ii) rationalize, and (iii) predict the effect of viral mutations on the binding affinity. Here, we investigate the allosteric signal triggered by the dissociation of the ACE2-RBD complex. To this end, we construct an Elastic Network Model (ENM), and we use the Structural Perturbation Method (SPM). Our key result is that complex dissociation opens the ACE2 substrate-binding cleft located away from the interface and that fluctuations of the ACE2 binding cleft are facilitated by RBD binding. These and other observations provide a structural and dynamical basis for the influence of SARS-CoV-2 on ACE2 enzymatic activity. In addition, we identify a conserved glycine (G502 in SARS-CoV-2) as a key participant in complex disassembly.

The Many Ways to Deal with STING.

The stimulator of interferon genes (STING) is an adaptor protein involved in the activation of IFN-ß and many other genes associated with the immune response activation in vertebrates. STING induction has gained attention from different angles such as the potential to trigger an early immune response against different signs of infection and cell damage, or to be used as an adjuvant in cancer immune treatments. Pharmacological control of aberrant STING activation can be used to mitigate the pathology of some autoimmune diseases. The STING structure has a well-defined ligand binding site that can harbor natural ligands such as specific purine cyclic di-nucleotides (CDN). In addition to a canonical stimulation by CDNs, other non-canonical stimuli have also been described, whose exact mechanism has not been well defined. Understanding the molecular insights underlying the activation of STING is important to realize the different angles that need to be considered when designing new STING-binding molecules as therapeutic drugs since STING acts as a versatile platform for immune modulators. This review analyzes the different determinants of STING regulation from the structural, molecular, and cell biology points of view.

Novel Polymyxin-Inspired Peptidomimetics Targeting the SARS-CoV-2 Spike:hACE2 Interface.

Though the bulk of the COVID-19 pandemic is behind, the search for effective and safe anti-SARS-CoV-2 drugs continues to be relevant. A highly pursued approach for antiviral drug development involves targeting the viral spike (S) protein of SARS-CoV-2 to prevent its attachment to the cellular receptor ACE2. Here, we exploited the core structure of polymyxin B, a naturally occurring antibiotic, to design and synthesize unprecedented peptidomimetics (PMs), intended to target contemporarily two defined, non-overlapping regions of the S receptor-binding domain (RBD). Monomers 1, 2, and 8, and heterodimers 7 and 10 bound to the S-RBD with micromolar affinity in cell-free surface plasmon resonance assays (KD ranging from 2.31 µM to 2.78 µM for dimers and 8.56 µM to 10.12 µM for monomers). Although the PMs were not able to fully protect cell cultures from infection with authentic live SARS-CoV-2, dimer 10 exerted a minimal but detectable inhibition of SARS-CoV-2 entry in U87.ACE2+ and A549.ACE2.TMPRSS2+ cells. These results validated a previous modeling study and provided the first proof-of-feasibility of using medium-sized heterodimeric PMs for targeting the S-RBD. Thus, heterodimers 7 and 10 may serve as a lead for the development of optimized compounds, which are structurally related to polymyxin, with improved S-RBD affinity and anti-SARS-CoV-2 potential.

Long Time Scale Ensemble Methods in Molecular Dynamics: Ligand-Protein Interactions and Allostery in SARS-CoV-2 Targets.

We subject a series of five protein-ligand systems which contain important SARS-CoV-2 targets, 3-chymotrypsin-like protease (3CLPro), papain-like protease, and adenosine ribose phosphatase, to long time scale and adaptive sampling molecular dynamics simulations. By performing ensembles of ten or twelve 10 µs simulations for each system, we accurately and reproducibly determine ligand binding sites, both crystallographically resolved and otherwise, thereby discovering binding sites that can be exploited for drug discovery. We also report robust, ensemble-based observation of conformational changes that occur at the main binding site of 3CLPro due to the presence of another ligand at an allosteric binding site explaining the underlying cascade of events responsible for its inhibitory effect. Using our simulations, we have discovered a novel allosteric mechanism of inhibition for a ligand known to bind only at the substrate binding site. Due to the chaotic nature of molecular dynamics trajectories, regardless of their temporal duration individual trajectories do not allow for accurate or reproducible elucidation of macroscopic expectation values. Unprecedentedly at this time scale, we compare the statistical distribution of protein-ligand contact frequencies for these ten/twelve 10 µs trajectories and find that over 90% of trajectories have significantly different contact frequency distributions. Furthermore, using a direct binding free energy calculation protocol, we determine the ligand binding free energies for each of the identified sites using long time scale simulations. The free energies differ by 0.77 to 7.26 kcal/mol across individual trajectories depending on the binding site and the system. We show that, although this is the standard way such quantities are currently reported at long time scale, individual simulations do not yield reliable free energies. Ensembles of independent trajectories are necessary to overcome the aleatoric uncertainty in order to obtain statistically meaningful and reproducible results. Finally, we compare the application of different free energy methods to these systems and discuss their advantages and disadvantages. Our findings here are generally applicable to all molecular dynamics based applications and not confined to the free energy methods used in this study.

Investigation on the binding behavior of human α1-acid glycoprotein with Janus Kinase inhibitor baricitinib: Multi-spectroscopic and molecular simulation methodologies.

Baricitinib is a Janus Kinase (JAK) inhibitor that is primarily used to treat moderately to severely active rheumatoid arthritis in adults and has recently been reported for the treatment of patients with severe COVID-19. This paper describes the investigation of the binding behavior of baricitinib to human α1-acid glycoprotein (HAG) employing a variety of spectroscopic techniques, molecular docking and dynamics simulations. Baricitinib can quench the fluorescence from amino acids in HAG through a mix of dynamic and static quenching, according to steady-state fluorescence and UV spectra observations, but it is mainly static quenching at low concentration. The binding constant (Kb) of baricitinib to HAG at 298 K was at the level of 104 M-1, indicating a moderate affinity of baricitinib to HAG. Hydrogen bonding and hydrophobic interactions conducted the main effect, according to thermodynamic characteristics, competition studies between ANS and sucrose, and molecular dynamics simulations. For the change in HAG conformation, the results of multiple spectra showed that baricitinib was able to alter the secondary structure of HAG as well as increase the polarity of the microenvironment around the Trp amino acid. Furthermore, the binding behavior of baricitinib to HAG was investigated by molecular docking and molecular dynamics simulations, which validated experimental results. Also explored is the influence of K+, Co2+, Ni2+, Ca2+, Fe3+, Zn2+, Mg2+ and Cu2+plasma on binding affinity.

Molecular Basis for the Selectivity of DHA and EPA in Sudlow's Drug Binding Sites in Human Serum Albumin with the Combined Use of NMR and Docking Calculations.

Medium- and long-chain saturated and unsaturated free fatty acids (FFAs) are known to bind to human serum albumin (HSA), the main plasma carrier protein. Atomic-level structural data regarding the binding mode in Sudlow's sites I (FA7) and II (FA4, FA3) of the polyunsaturated ω-3 fatty acids docosahexaenoic acid (DHA) and eicosapentaenoic acid (EPA), however, are largely unknown. Herein, we report the combined use of saturation transfer difference (STD) and Interligand NOEs for Pharmacophore Mapping (INPHARMA) NMR techniques and molecular docking calculations to investigate the binding mode of DHA and EPA in Sudlow's sites Ι and ΙΙ of HSA. The docking calculations and the significant number of interligand NOEs between DHA and EPA and the drugs warfarin and ibuprofen, which are stereotypical ligands for Sudlow's sites I and II, respectively, were interpreted in terms of competitive binding modes and the presence of two orientations of DHA and EPA at the binding sites FA7 and FA4. The exceptional flexibility of the long-chain DHA and EPA and the formation of strongly folded structural motives are the key properties of HSA-PUFA complexes.

Computational and Enzymatic Studies of Sartans in SARS-CoV-2 Spike RBD-ACE2 Binding: The Role of Tetrazole and Perspectives as Antihypertensive and COVID-19 Therapeutics.

This study is an extension of current research into a novel class of synthetic antihypertensive drugs referred to as "bisartans", which are bis-alkylated imidazole derivatives bearing two symmetric anionic biphenyltetrazoles. Research to date indicates that bisartans are superior to commercially available hypertension drugs, since the former undergo stronger docking to angiotensin-converting enzyme 2 (ACE2). ACE2 is the key receptor involved in SARS-CoV-2 entry, thus initiating COVID-19 infection and in regulating levels of vasoactive peptides such as angiotensin II and beneficial heptapeptides A(1-7) and Alamandine in the renin-angiotensin system (RAS). In previous studies using in vivo rabbit-iliac arterial models, we showed that Na+ or K+ salts of selected Bisartans initiate a potent dose-response inhibition of vasoconstriction. Furthermore, computational studies revealed that bisartans undergo stable binding to the vital interfacial region between ACE2 and the SARS-CoV-2 "receptor binding domain" (i.e., the viral RBD). Thus, bisartan homologs are expected to interfere with SARS-CoV-2 infection and/or suppress disease expression in humans. The primary goal of this study was to investigate the role of tetrazole in binding and the network of amino acids of SARS-CoV-2 Spike RBD-ACE2 complex involved in interactions with sartans. This study would, furthermore, allow the expansion of the synthetic space to create a diverse suite of new bisartans in conjunction with detailed computational and in vitro antiviral studies. A critical role for tetrazole was uncovered in this study, shedding light on the vital importance of this group in the binding of sartans and bisartans to the ACE2/Spike complex. The in silico data predicting an interaction of tetrazole-containing sartans with ACE2 were experimentally validated by the results of surface plasmon resonance (SPR) analyses performed with a recombinant human ACE2 protein.

LigBind: Identifying Binding Residues for Over 1000 Ligands with Relation-Aware Graph Neural Networks.

Identifying the interactions between proteins and ligands is significant for drug discovery and design. Considering the diverse binding patterns of ligands, the ligand-specific methods are trained per ligand to predict binding residues. However, most of the existing ligand-specific methods ignore shared binding preferences among various ligands and generally only cover a limited number of ligands with a sufficient number of known binding proteins. In this study, we propose a relation-aware framework LigBind with graph-level pre-training to enhance the ligand-specific binding residue predictions for 1159 ligands, which can effectively cover the ligands with a few known binding proteins. LigBind first pre-trains a graph neural network-based feature extractor for ligand-residue pairs and relation-aware classifiers for similar ligands. Then, LigBind is fine-tuned with ligand-specific binding data, where a domain adaptive neural network is designed to automatically leverage the diversity and similarity of various ligand-binding patterns for accurate binding residue prediction. We construct ligand-specific benchmark datasets of 1159 ligands and 16 unseen ligands, which are used to evaluate the effectiveness of LigBind. The results demonstrate the LigBind's efficacy on large-scale ligand-specific benchmark datasets, and it generalizes well to unseen ligands. LigBind also enables accurate identification of the ligand-binding residues in the main protease, papain-like protease and the RNA-dependent RNA polymerase of SARS-CoV-2. The web server and source codes of LigBind are available at http://www.csbio.sjtu.edu.cn/bioinf/LigBind/ and https://github.com/YYingXia/LigBind/ for academic use.

Analysis of the ARTIC V4 and V4.1 SARS-CoV-2 primers and their impact on the detection of Omicron BA.1 and BA.2 lineage-defining mutations.

The ARTIC protocol uses a multiplexed PCR approach with two primer pools tiling the entire SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2) genome. Primer pool updates are necessary for accurate amplicon sequencing of evolving SARS-CoV-2 variants with novel mutations. The suitability of the ARTIC V4 and updated V4.1 primer scheme was assessed using whole genome sequencing of Omicron from clinical samples using Oxford Nanopore Technology. Analysis of Omicron BA.1 genomes revealed that 93.22 % of clinical samples generated improved genome coverage at 50× read depth with V4.1 primers when compared to V4 primers. Additionally, the V4.1 primers improved coverage of BA.1 across amplicons 76 and 88, which resulted in the detection of the variant-defining mutations G22898A, A26530G and C26577G. The Omicron BA.2 sub-variant (VUI-22JAN-01) replaced BA.1 as the dominant variant by March 2022, and analysis of 168 clinical samples showed reduced coverage across amplicons 15 and 75. Upon further interrogation of primer binding sites, a mutation at C4321T [present in 163/168 (97 %) of samples] was identified as a possible cause of complete dropout of amplicon 15. Furthermore, two mutations were identified within the primer binding regions for amplicon 75: A22786C (present in 90 % of samples) and C22792T (present in 12.5 % of samples). Together, these mutations may result in reduced coverage of amplicon 75, and further primer updates would allow the identification of the two BA.2-defining mutations present in amplicon 75: A22688G and T22679C. This work highlights the need for ongoing surveillance of primer matches as circulating variants evolve and change.

A universal fluorescence polarization high throughput screening assay to target the SAM-binding sites of SARS-CoV-2 and other viral methyltransferases.

SARS-CoV-2 has caused a global pandemic with significant humanity and economic loss since 2020. Currently, only limited options are available to treat SARS-CoV-2 infections for vulnerable populations. In this study, we report a universal fluorescence polarization (FP)-based high throughput screening (HTS) assay for SAM-dependent viral methyltransferases (MTases), using a fluorescent SAM-analogue, FL-NAH. We performed the assay against a reference MTase, NSP14, an essential enzyme for SARS-CoV-2 to methylate the N7 position of viral 5'-RNA guanine cap. The assay is universal and suitable for any SAM-dependent viral MTases such as the SARS-CoV-2 NSP16/NSP10 MTase complex and the NS5 MTase of Zika virus (ZIKV). Pilot screening demonstrated that the HTS assay was very robust and identified two candidate inhibitors, NSC 111552 and 288387. The two compounds inhibited the FL-NAH binding to the NSP14 MTase with low micromolar IC50. We used three functional MTase assays to unambiguously verified the inhibitory potency of these molecules for the NSP14 N7-MTase function. Binding studies indicated that these molecules are bound directly to the NSP14 MTase with similar low micromolar affinity. Moreover, we further demonstrated that these molecules significantly inhibited the SARS-CoV-2 replication in cell-based assays at concentrations not causing cytotoxicity. Furthermore, NSC111552 significantly synergized with known SARS-CoV-2 drugs including nirmatrelvir and remdesivir. Finally, docking suggested that these molecules bind specifically to the SAM-binding site on the NSP14 MTase. Overall, these molecules represent novel and promising candidates to further develop broad-spectrum inhibitors for the management of viral infections.

De novo design of protein interactions with learned surface fingerprints.

Physical interactions between proteins are essential for most biological processes governing life1. However, the molecular determinants of such interactions have been challenging to understand, even as genomic, proteomic and structural data increase. This knowledge gap has been a major obstacle for the comprehensive understanding of cellular protein-protein interaction networks and for the de novo design of protein binders that are crucial for synthetic biology and translational applications2-9. Here we use a geometric deep-learning framework operating on protein surfaces that generates fingerprints to describe geometric and chemical features that are critical to drive protein-protein interactions10. We hypothesized that these fingerprints capture the key aspects of molecular recognition that represent a new paradigm in the computational design of novel protein interactions. As a proof of principle, we computationally designed several de novo protein binders to engage four protein targets: SARS-CoV-2 spike, PD-1, PD-L1 and CTLA-4. Several designs were experimentally optimized, whereas others were generated purely in silico, reaching nanomolar affinity with structural and mutational characterization showing highly accurate predictions. Overall, our surface-centric approach captures the physical and chemical determinants of molecular recognition, enabling an approach for the de novo design of protein interactions and, more broadly, of artificial proteins with function.

SARS-CoV-2 viral spike G614 mutation exhibits higher case fatality rate.

AIM: The COVID-19 pandemic is caused by infection with the SARS-CoV-2 virus. The major mutation detected to date in the SARS-CoV-2 viral envelope spike protein, which is responsible for virus attachment to the host and is also the main target for host antibodies, is a mutation of an aspartate (D) at position 614 found frequently in Chinese strains to a glycine (G). We sought to infer health impact of this mutation. RESULT: Increased case fatality rate correlated strongly with the proportion of viruses bearing G614 on a country by country basis. The amino acid at position 614 occurs at an internal protein interface of the viral spike, and the presence of G at this position was calculated to destabilise a specific conformation of the viral spike, within which the key host receptor binding site is more accessible. CONCLUSION: These results imply that G614 is a more pathogenic strain of SARS-CoV-2, which may influence vaccine design. The prevalence of this form of the virus should also be included in epidemiologic models predicting the COVID-19 health burden and fatality over time in specific regions. Physicians should be aware of this characteristic of the virus to anticipate the clinical course of infection.

Effect of adjuvanting RBD-dimer-based subunit COVID-19 vaccines with Sepivac SWE™.

Protein subunit vaccines have been widely used to combat infectious diseases, including the current COVID-19 pandemic. Adjuvants play the key role in shaping the quality and magnitude of the immune response to protein and inactivated vaccines. We previously developed a protein subunit COVID-19 vaccine, termed ZF2001, based on an aluminium hydroxide-adjuvanted tandem-repeat dimeric receptor-binding domain (RBD) of the viral spike (S) protein. Here, we described the use of a squalene-based oil-in-water adjuvant, Sepivac SWE™ (abbreviated to SWE), to further improve the immunogenicity of this RBD-dimer-based subunit vaccines. Compared with ZF2001, SWE adjuvant enhanced the antibody and CD4+ T-cell responses in mice with at least 10 fold of dose sparing compared with ZF2001 adjuvanted with aluminium hydroxide. SWE-adjuvanted vaccine protected mice against SARS-CoV-2 challenge. To ensure adequate protection against the currently circulating Omicron variant, we evaluated this adjuvant in combination with Delta-Omicron chimeric RBD-dimer. SWE significantly increased antibody responses compared with aluminium hydroxide adjuvant and afforded greater neutralization breadth. These data highlight the advantage of emulsion-based adjuvants to elevate the protective immune response of protein subunit COVID-19 vaccines.

Identification of Four New Chemical Series of Small Drug-Like Natural Products as Potential Neuropilin-1 Inhibitors by Structure-Based Virtual Screening: Pharmacophore-Based Molecular Docking and Dynamics Simulation.

Neuropilin-1 (NRP-1), a surface transmembrane glycoprotein, is one of the most important co-receptors of VEGF-A165 (vascular endothelial growth factor) responsible for pathological angiogenesis. In general, NRP-1 overexpression in cancer correlates with poor prognosis and more tumor aggressiveness. NRP-1 role in cancer has been mainly explained by mediating VEGF-A165-induced effects on tumor angiogenesis. NRP-1 was recently identified as a co-receptor and an independent gateway for SARS-CoV-2 through binding subunit S2 of Spike protein in the same way as VEGF-A165. Thus, NRP-1 is of particular value as a target for cancer therapy and other angiogenesis-dependent diseases as well as for SARS-CoV-2 antiviral intervention. Herein, The Super Natural II, the largest available database of natural products (∼0.33 M), pre-filtered with drug-likeness criteria (absorption, distribution, metabolism and excretion/toxicity), was screened against NRP-1. NRP-1/VEGF-A165 interaction is one of protein-protein interfaces (PPIs) known to be challenging when approached in-silico. Thus, a PPI-suited multi-step virtual screening protocol, incorporating a derived pharmacophore with molecular docking and followed by MD (molecular dynamics) simulation, was designed. Two stages of pharmacophorically constrained molecular docking (standard and extra precisions), a mixed Torsional/Low-mode conformational search and MM-GBSA ΔG binding affinities calculation, resulted in the selection of 100 hits. These 100 hits were subjected to 20 ns MD simulation, that was extended to 100 ns for top hits (20) and followed by post-dynamics analysis (atomic ligand-protein contacts, RMSD, RMSF, MM-GBSA ΔG, Rg, SASA and H-bonds). Post-MD analysis showed that 19 small drug-like nonpeptide natural molecules, grouped in four chemical scaffolds (purine, thiazole, tetrahydropyrimidine and dihydroxyphenyl), well verified the derived pharmacophore and formed stable and compact complexes with NRP-1. The discovered molecules are promising and can serve as a base for further development of new NRP-1 inhibitors.

Probing Ligand Binding Sites on Large Proteins by Nuclear Magnetic Resonance Spectroscopy of Genetically Encoded Non-Canonical Amino Acids.

N6-(((trimethylsilyl)-methoxy)carbonyl)-l-lysine (TMSK) and N6-trifluoroacetyl-l-lysine (TFAK) are non-canonical amino acids, which can be installed in proteins by genetic encoding. In addition, we describe a new aminoacyl-tRNA synthetase specific for N6-(((trimethylsilyl)methyl)-carbamoyl)-l-lysine (TMSNK), which is chemically more stable than TMSK. Using the dimeric SARS-CoV-2 main protease (Mpro) as a model system with three different ligands, we show that the 1H and 19F nuclei of the solvent-exposed trimethylsilyl and CF3 groups produce intense signals in the nuclear magnetic resonance (NMR) spectrum. Their response to active-site ligands differed significantly when positioned near rather than far from the active site. Conversely, the NMR probes failed to confirm the previously reported binding site of the ligand pelitinib, which was found to enhance the activity of Mpro by promoting the formation of the enzymatically active dimer. In summary, the amino acids TMSK, TMSNK, and TFAK open an attractive path for site-specific NMR analysis of ligand binding to large proteins of limited stability and at low concentrations.

Increased flexibility of the SARS-CoV-2 RNA-binding site causes resistance to remdesivir.

Mutations continue to accumulate within the SARS-CoV-2 genome, and the ongoing epidemic has shown no signs of ending. It is critical to predict problematic mutations that may arise in clinical environments and assess their properties in advance to quickly implement countermeasures against future variant infections. In this study, we identified mutations resistant to remdesivir, which is widely administered to SARS-CoV-2-infected patients, and discuss the cause of resistance. First, we simultaneously constructed eight recombinant viruses carrying the mutations detected in in vitro serial passages of SARS-CoV-2 in the presence of remdesivir. We confirmed that all the mutant viruses didn't gain the virus production efficiency without remdesivir treatment. Time course analyses of cellular virus infections showed significantly higher infectious titers and infection rates in mutant viruses than wild type virus under treatment with remdesivir. Next, we developed a mathematical model in consideration of the changing dynamic of cells infected with mutant viruses with distinct propagation properties and defined that mutations detected in in vitro passages canceled the antiviral activities of remdesivir without raising virus production capacity. Finally, molecular dynamics simulations of the NSP12 protein of SARS-CoV-2 revealed that the molecular vibration around the RNA-binding site was increased by the introduction of mutations on NSP12. Taken together, we identified multiple mutations that affected the flexibility of the RNA binding site and decreased the antiviral activity of remdesivir. Our new insights will contribute to developing further antiviral measures against SARS-CoV-2 infection.

Interactions of angiotensin-converting enzyme-2 (ACE2) and SARS-CoV-2 spike receptor-binding domain (RBD): a structural perspective.

BACKGROUND: Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) has caused millions of infections and deaths worldwide since its discovery in late 2019 in Wuhan, China. The receptor-binding domain (RBD) of the SARS-CoV-2 spike protein binds to the human angiotensin-converting enzyme-2 (ACE2) receptor, a critical component of the renin-angiotensin system (RAS) that initiates the viral transmission. Most of the critical mutations found in SARS-CoV-2 are associated with the RBD of the spike protein. These mutations have the potential to reduce the efficacy of vaccines and neutralizing antibodies. METHODS: In this review, the structural details of ACE2, RBD and their interactions are discussed. In addition, some critical mutations of RBD and their impact on ACE2-RBD interactions are also discussed. CONCLUSION: Preventing the interaction between Spike RBD and ACE2 is considered a viable therapeutic strategy since ACE2 binding by RBD is the first step in virus infection. Because the interactions between the two entities are critical for both viral transmission and therapeutic development, it is essential to understand their interactions in detail.